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ABSTRACT 

For more than 15 years, since the days of the Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton Gamma-Ray 
Observatory (CGRO; 1991-2000), it has remained an open question why the prominent blazar 3C345 was not reliably detected at y-ray ener- 
gies >20MeV. Recently a bright y-ray source (OFGL J1641.4+3939/1FGL J1642.5+3947), potentially associated with 3C345, was detected by 
the Large Area Telescope (LAT) on Fermi. Multiwavelength observations from radio bands to X-rays (mainly GASP-WEBT and Swift) of possible 
counterparts (3C345, NRA0 512, B3 1640+396) were combined with 20 months of FerraZ-LAT monitoring data (August 2008 - April 2010) to 
associate and identify the dominating y-ray emitting counterpart of 1FGLJ1642. 5+3947. The source 3C 345 is identified as the main contributor 
for this y-ray emitting region. However, after November 2009 (15 months), a significant excess of photons from the nearby quasar NRA0512 
started to contribute and thereafter was detected with increasing y-ray activity, possibly adding flux to 1FGLJ1642. 5+3947. For the same time 
period and during the summer of 2010, an increase of radio, optical and X-ray activity of NRA0512 was observed. No y-ray emission from 
B3 1640+396 was detected. 
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The quasar 3C 345 is known as a prominent variable soi rce from 
radio- to X-ray bands (< lOkeV). It is particularly bright at ra- 
dio wavelengths and has an extended radio structure diat is ob- 
servable from sub-pc to kpc scales which is archety p cal for a 
relati vistic blazar j et (e.g. 'Lobanov & Zensus ni999l: IMlgaard 
et al. Il989l). Its jet has an apparent opening angle of ~ 13° 



(Pushkarev et al. 2009). IUnwin et alJ (11994) concluded that the 
observed high-energy emission of this source is caus(;d by in- 
verse Compton up-scattering of lower energy photons by rela- 
tivistic electrons of the pc-scale radio jet. However, it i^emained 
a puzzle why 3C 345 was not detected at energies above 20 MeV, 
while 3C279 with similar apparent properties is a plrominent 
source at y-ray energies. | 

The high-redshift quasar NRAO 512, 0.5° west of 3C 345, is 
known as a compact radio source with a two-sided jet morphol- 
ogyon kpc scales and with a flat radio spectrum (Kovalev et al. 



Il999l) . The observed VLBI radio structure revealed a 



wide ap- 



parent jet opening angle of ~ 50° jPushkarev et al.ll200 3), which 



implies a small viewing angle to the line of sight, with its emis- 
sion directly beamed at us. 

The Energetic Gamm a-Ray Experiment Telescope (EGRET; 
[Thompson et al.l Il993h on board the Compton Gamma Ray 
Observatory (CGRO) spacecraft did not observe a significant 
excess of y-ray photons consistent with 3C 345 (3EG; iHartman 



et al.| [T999h . covering the first four phases of the CGRO mis- 
sion [1991-1995). This lack of detectable flux of y-rays above 
20 MeV from 3C 345 was attributed to the high detection 
threshold of EGRET. A source consistent with the position of 
3C 345, labeled as EGR J 1642+3940. was found bv Casandiian 
& Grenier (2008). This source was seen almost entirely in a sin- 
gle two-week viewing period with an excess of ~5.8cr, from 
April 23rd - May 7th, 1996 (vp5190), which was after the 
data used for the 3EG. Unfortunately, the angular resolution 
of E<jRET was too poor to distinguish contributions from the 
other two nearby candidates NRAO 5 12 and B3 1640+396, see 
Table; [1] Between October 1st and October 4th, 1993, a simi- 
lar source was marginally detected by EGRET (~2.1cr), but its 
centrpid was lo calized closer t o Mrk 501. In exactly the same 
view ng period. iKataoka et al.l (Il999h claimed EGRET detected 
a flar; from Mrk501, which lies 2.5° east of 3C 345. This would 
then be the only time that either of those sources was seen by 
EGRET - a strange coincidence that led Casandjian & Grenier to 
note that the EGR J 1642+3940 source association with 3C345 
is "not clear". 



On. 

lOljt I 
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I June 1 1th, 2008, the Fermi Gamma-ray Space Telescope 
(shoijt Fermi) was launched, equipped with the Large Area 
Telescope (LAT) instrument (iAtwood et alJ2009 ) providing con- 
tinuous all-sky monitoring of y-ray emission as a default mode 
of op bration. The LAT, compared with EGRET, has an improved 
sensitivity by a factor of ~ 25 and an improved point-source 
localization. After 20 months of LAT operation within the re- 
gion of interest, a localization improvement (comparing radii 
of 95% error circles) of a factor of ~ 14 was reached compared 
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Table 1. Radio positions and redshifts of candidate y-ray coun- 
terparts. 



Name(s) 


3C345 


NRA0 512 


B3 1640+396 
CLASSJ1641+3935 


R.A. 


16H2'"5S810 


16M0™29';633 


16H1"^47^:540 


Dec. 


+39°48'36'.'99 


+39°46'46'.'03 


+39°35'03'.'33 


z 


0.5928 


1.666 


0.540 


Refs. 


(1,3) 


(1,4) 


(2,5) 



Re fs. - Tli e numerals in tlie references row translate to: 

(l)| Fevetal.(2004);(2 ) M yers et al. (2003); ( 3) Marziani et alj ( 119961) ; 
(4) 



Stickeletal.l ( ll989h : (5) lBrinkmann et alJ ( ll995l) 



Notes: Name(s) - Common source names. 

R.A./Dec. - Right ascension/declination (J2000) of the radio source 

localization. 

z - Source redshift. 



to EGRET. These improvements provided a new opportunity to 
probe 3C 345 for y-ray emission above 100 Me V. 

The three-month bright Fermi source list (OFGL) for the 
dataset collected between August and October 2008 reported 
a y-ray source (0FGLJ1641. 4+3939) near the position of 
3C345. It was then associated with the quasar B3 1640+39 6 
(CLASS J1641+3935) at low confidence (lAbdo et al.ll2009bh . 
The bright OFGL source was listed in the first Fermi-LK^ so urce 
catalog (IFGL) as IFGL J1642.5+3947 (lAbdoetall2010al) . fur- 
ther constraining its localization. An association with 3C345 
was still not possible with high confidence based on Fermi-LP^ 
data alone, nevertheless, IFGL J1642. 5+3947 was listed in the 
catalog as being afli fiated with its most likely counterpart 3C 345 
(lAbdo et al.ll2010bh . 

In October 2009 a GeV flar e was detected from th e 
vicinity of OFGL J 164 1.4+3939 (iReves & Cheund l2009h . 
Simultaneously, the GASP-WEBT collaboratio n reported on in- 
crease d optical and radio activity of 3C 345 (iLarionov et al.l 
l2009l) . This event triggered additional observations at different 
frequencies, which made it possible to establish an identification 
through coiTelated multiwavelength activity. 

Here we present 20 months of monitoring observations 
by Fermi LAT of the region around 3C345, NRA0512 and 
B3 1640+396, as well as radio, optical, UV, and X-ray obser- 
vations of these sources. Section 2 discusses the observations 
and data reduction methods applied to obtain calibrated datasets 
for the subsequent analysis. Section 3 discusses the results of 
Ferm/-LAT observations and identification through multiwave- 
length counterparts of events in the light curve. Section 4 dis- 
cusses and summarizes the findings of this paper. 

2. Observations and data reduction metliods 

2.1. y-ray monitoring 

The Fermi LAT (lAtwood et al.ll2009h is a pair-conversion tele- 
scope designed to cover the energy range from 20 MeV to greater 
than 300 GeV. The IFGL catalog, based on 11 months of mon- 
itoring data, lists the y-ray source 1FGLJ1642. 5+3947 (R.A. 
16''42"'33.55 \ Dec. +39°47'38 .8", 95% position error radius 
r95=0.056°: .Abdo et al.ll2010ah . which lies in close vicinity of 
the three quasars 3C 345, NRAO 512 and B3 1640+396, listed in 
Table [U 

The Fermi-LAiT data on the field in which 
IFGL J1642. 5+3947 is located were processed and ana- 
lyzed using the Fermi Science Tools software package (version: 



v9rl5p2Pl: lAbdo et ani2010ah . The photons were extracted from 
a region of interest (ROI) centered on the radio position of 
3C 345 within a radius of 15° and a time range of 20 months, 
August 4th, 2008, to April 4th, 2010 (in MElQ 239557417 - 
292088627 s). The general data processing procedur es applied 
were similar to those described in Abdo et aTl (l2009bl) and Abdo 



et al. (I2010ah . 

To ensure that the collected events considered in this analy- 
sis have a high probability of being photons, the "diffuse class" 
selection was applied. Furthermore, events above the zenith an- 
gle 105° were removed to avoid a significant contamination by 
y rays p roduced from cosmic -ray interactions in Earth's atmo- 
spher^ (lAbdo et al.ll2009c) . The P6_V3_DIFFUSE instrument 
response functions (IRF) were used (see lRando et al.ll2009l) . 

The remaining y-ray-like events wer e analyzed using a 
maximum-likelihood approach dMattox et al. 1996) to localize 
y-ray sources and extract their spectra. As in the analysis of 
IFGL, two sets of standard LAT analysis tools were used, and 
their results compared for consistency. The first set, based on 
gtfindsrc and gtlike tools, employs an unbinned likelihood 
analysis in which the expected response of the LAT is evaluated 
separately for each photon, given its energy and direction with 
respect to the detector axes. The second toolset, pointlike, 
uses a simpler, binned analysis methodology, where photons are 
binned spatially and in energy, and an averaged set of response 
functions is applied. The two packages give very similar results, 
with the pointlike package running much more quickly. 

In both likelihood methodologies the y-ray emission in the 
ROI is assumed to arise from two diffuse components and a num- 
ber of point sources. Each point source is modeled spatially as a 
delta function and spectrally as a power law. The statistical sig- 
nificance of a potential source is estimated by forming a "test 
statistic" (TS), 



TS =2(lnLi -InLo), 



(1) 



where Li is the value of the likelihood obtained in the "source" 
hypothesis, with the source present in the model of the ROI and 
Lq is the value in the alternate, "null" hypothesis, that the source 
is not real (and hence not included in the ROI model). If the 
null hypothesis is true, and assumi ng certai n applicability condi- 
tions are met, Wilks' theorem ( Wilks 1938) states that TS is dis- 
tributed as x^ (N), where A^ is the number of additional parame- 
ters optimized in the source model leading to Li . Because each 
point source has four free parameters (two for its position and 
two for its spectrum), we expect that TS ~ x^ (4), from which 
its significance can be evaluated. A source-detection threshold 
of TS=25 was adopted in IFGL, which corresponds to a false 
detection probability of 5 x 10"^ if TS ~ x^ (4). 

The gtlike model used for localization in this paper in- 
cludes the IFGL J1642. 5+3947 point-sourc e component and al l 
other point sources from the IFGL catalog dAbdo et al.ll2010a) 
within a 20° radius around the radio position of 3C 345. Further 
analysis was performed using two point-source components in- 
stead of one to fit the observed data. The model parameters for 
sources outside the 15° ROI were fixed to their respective IFGL 
values. The model included 20 or 21 point sources with free 
parameters (flux and spectral index), depending on whether a 



' http://fermi.gsfc. nasa.gov/.. ./documentation/Cicerone 
- Mission Elapsed Time (MET) is measured in seconds from 2001.0 
UT, including leap seconds. 

■^ The Earth's limb lies at a zenith angle of ~ 11 3° at the 565 km nearly 
circular orbit of Fermi. The zenith angle is defined as the angle of a 
photon's apparent origin to the Earth-spacecraft vector. 
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single point source was fitted for 1FGLJ1642. 5+3947 (20 point 
sources) or whether two source components were added to take 
into account the separate contributions of 3C 345 and NRAO 512 
(21 point sources). The background component of Galactic 
diffuse emission (gll_iem_v02 .fit) and an isotropic compo- 
nent (isotropic_iem_v02.txt), both of which are the stan- 
dard models available from the Fermi Science Support Center 
(FSSCJj, were added as well. The isotropic component includes 
both the contribution from the extragalactic diffuse emission and 
from the residual, charged particle backgrounds. 

To produce light curves, the standard tool gtlike was used, 
applying an unbinned spectral likelihood analysis as described 
above. A gtlike input source model was constructed based on 
the improved 20-month spectral fits obtained from the localiza- 
tion analysis. Only sources that were detected with a TS > 25 in 
the 20-month integration were included in the final gtlike in- 
put source model to construct a light curve. The data were split 
into regular time intervals, each integrating over periods of 2, 5 
and 30 days. Then an unbinned spectral likelihood analysis was 
performed on each integrated dataset with photon indices fixed 
to the 20-month average values. For time periods with a signifi- 
cance of the detection of TS < 5 on 5- and 30-day scales, 2cr up- 
per limits wer e calculated following the likelihood profile (Rolke 
et al. 2005), where the upper limit is determined by increasing 
the flux from the maximum likelihood value until the log likeli- 
hood decreases by 2.0 (for the part icular implementation in case 
of LAT data see lAbdo et al.ll2010ah . All light curves make use of 
the collected data in the energy range of 0. l-300GeV. Different 
time binnings were chosen to study different aspects of the time 
evolution of the source flux. To study the longer term behavior, 
30-day bins were used (Section |3]2TTJ, monthly bins being rela- 
tively standard for Fermi analysis. To assess faster flaring behav- 
ior, and in particular to estimate when the peak flaring occurs, a 
two-day time binning was selected (Section [3.2.2b . Finally, the 
best compromise between short and long time bins is the light 
curve with five-day integration periods (Section [33T l. 

2.2. X-ray observations 

Triggered by t he reported y-ray flare of OFGL J 164 1.4-1-3939 in 
October 2009 dReves & Cheun^l2009t) . X-ray observations of all 
three candidate counterpart sour ces were performed by the X- 
ray Telescope (XR T, 0.2-10keV ;'Burrow s et alJl2005h on board 
the Swift satellite dGehrels et ah.2004,) . The Ust of Swift-XRT 
observations together with the spectral fit results is reported in 
Tabled 

The XRT data were processed with standard procedures 
(xrtpipeline vO.12.4) to calibrate the observations. The filter- 
ing and screening criteria were applied by means of the FTOOLS 
in the Heasoft package version 6.qj. Given the low fluxes of the 
sources during observations (< 0.5 counts s"' in the 0.2-lOkeV 
range), only photon counting (PC) data were considered for our 
analysis with XR T grades 0-12 sele cted (according to Swift 
nomenclature, see lBurrows et alj|2005 !). The source events were 
extracted in a circle of radius 15-20 pixels (covering 85-90% 
of the XRT point spread function, PSF, at 1 .5 keV) around the 
source, depending on its intensity. The background was esti- 
mated from a nearby source-free circular region of 40-pixel ra- 



"* http://fermi. gsfc.nasa.gov/.. ./Iat/BackgroundModels.html| 
"■ xrtpipeline and FTOOLS are part of the Heasoft software 
package: http://heasarc.gsfc.nasa.gov/Ilieasoft/, a detailed discussion 
on XRT data analysis can be found in the XRT User's Guide: 
[http: //heasarc. g sfc.nasa.gov/.. ./xrt-swguide-vl-l.pdf II 



dius. Ancillary response files were generated with xrtmkarf to 
account for different extraction regions, vignetting and PSF cor- 
rections. The spectral redistribution matrices vOll in the calibra- 
tion database maintained by HEASARC were used. 

The adopted energy range for spectral fitting is 0.3-lOkeV. 
All data were rebinned with a minimum of 20 counts per energy 
bin to use the x^ minimization fitting technique, except when the 
statistics were p oor (< 200 counts), in which case the Cash statis- 
tic (Cash 1 19791) was used and the data were binned to 1 counts 
bin '. The Cash statistic is based on a likelihood ratio test and 
is widely used for parameter estimation in photon counting ex- 
periments. Swift-XSH uncertainties are given at 90% confidence 
level for each parameter, unless stated otherwise. 

Spectral fitting was performed using XSPEC version 12.5.1 
with an absorbed power law model. For the Galactic absorption 
a fixed HI column density ne of 1.14 ■ 10^" cm"^ was used for 
3C 345, 1.17-102"cm-2 for B3 1640-H396 and 1.07-1020 cm-2 for 
NRAO 512 , obtained from the L eiden/ Argentine/Bonn Galactic 
H I survey (iKalberla et al.ll2005l) . 

2.3. UV, optical and mm observations 

In the case of 3C 345 a dense optical monitoring was per- 
formed by the GLAST-AGILE Support Program (GASP) of 
the Whole Earth Blazar Telescope (WEBT) collaboration. The 
GASP-WEBT is performing long-term monitoring of 28 y-ray 
loud blazars in the o ptical , near-infrared, mm and radio bands 
dVillata et alj l2008l l2009h . The GASP has been closely fol- 
lowing a strong opti cal outburst in 3C 345 reported during the 
second half of 2009 dLarionov et al.l l2009). The optical GASP 
data for this paper were acquired at the following observatories: 
Abastumani (Georgia), Calar Alto (Spain, within the MAPCAT 
program), Crimean (Ukraine), Lowell (Perkins, USA), Lu-lin 
(SLT, Taiwan), Sabadell (Spain), St. Petersburg (Russia), and 
Tijarafe (Spain). 

Magnitude calibration was performed with respect to a com- 
mon choice of reference stars in t he field of the source obtaine d 
from the photometric sequence bv iGonzal ez-Perez et alj d2001h . 
In addition to the GASP monitoring data, optical observations 
performed by the Special Astrophysical Observatory (SAO) of 
the Russian Academy of Sciences, using the SAO Zeiss 1 -m tele- 
scope in BVR bands were included. The observational strategy 
and calibration procedures applied were described in detail by 
Bvchkova et al (2004) . In Section [373] the R band magnitudes 
were converted to the linear flux density scale to compare light 
curves between wavelengths. The magnitudes were corrected for 
a Galactic extinction of 0.036 mag, for t he magnitude to flux 
conversion the zero magnitude fluxes of JBessell et alJ dl998h 
were used. 

Two R-band observations of 3C 345 and B3 1640H-396 were 
conducted on January 26th, 2009 (14x180 s exposure), and on 
February 13th, 2009 (13x180 s exposure), using a robotic 200- 
mm telescope of the Tzec Maun observatory (Mayhill, NM, 
USA). Magnitude calibration was performed in a manner sim- 
ilar to the GASP observations. 

Observations of 3C 345 and NRAO 512 were performed with 
the Ultravi olet/Optical Telesco pe (UVOT) on board the Swift 
spacecraft dRoming et alJl20()5h . Optical p hotometry of 3C345 
was calibrated using comparison stars from [Smith et alJ dl985h . 
Calibrat ion of the U V data was performed following the method 
of P oole et al.[ d2008l\ Comparison stars in the vicinity of 3C 345 
could not be used to calibrate optical photometry of NRAO 512 
with Swift UVOT because of its narrow field of view. Instead, the 
magnitude scale was calibrated using Sloan Digital Sky Survey 
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Table 2. Swift-XRT observations after the reported flare in October 2009 of the three candidate sources (3C345, NRA0 512 and 
B3 1640+396). 



Source 


Obs. Date (*) 


Exp. 


Flux 


Index 


3C 345: 


2009-10-04(5109) 


4.2 


6.8 ± 0.5 


1.68 + 0.11 




2009-10-06" (5111) 


1.9 


4.7 ± 0.7 


1.66 + 0.17 




2009-10-08(5113) 


2.1 


7.1 ±0.8 


1.75 + 0.16 




2009-10-09(5114) 


2.2 


6.4 ± 0.7 


1.62 + 0.15 




2009-10-18 (5123) 


3.9 


6.2 ± 0.5 


1.67 + 0.12 




2009-11-01 (5137) 


4.1 


5.5 ± 0.5 


1.77 + 0.11 




2009-11-15 (5151) 


4.0 


5.9 ± 0.5 


1.88 + 0.11 




2009-11-29(5165) 


4.1 


5.7 ± 0.5 


1.71 +0.12 




2010-03-06/07 (5262/3) 


2.6 


5.6 + 0.6 


1.84 + 0.15 




2010-03-09" (5265) 


1.4 


5.6 + 0.8 


1.78 + 0.18 




2010-04-09" (5296) 


2.0 


5.5 + 0.7 


1.88 + 0.16 




2010-08-18 (5427) 


4.3 


5.4 + 0.5 


1.54 + 0.11 


NRA0 512: 


2009-10-08" (5113) 


2.1 


0.6 !»i 


1.5 + 0.5 




2009-10-15" (5120) 


2.1 


0.5 + 0.2 


2.2 + 0.6 




2010-08-06/07" (5415/6) 


5.5 


0.7 + 0.1 


2.0 + 0.3 


B3 1640+396: 


2009-12-01" (5167) 


1.1 


7 *"-* 


1.8 + 0.7 



Notes: Obs. Date (*) - Date of Swift-XKI observation (*=JD-2450000). 

Exp. - Swift-XKI effective exposure time in ks. 

Flux - Integrated unabsorbed photon flux 0.3-10 keV in units of 10"'- erg cm"- s" 

" Cash statistic used for observations with less than 200 photons. 



(SDSS:lAbazaiian et alJ l2069h stars in the field and transforma- 
tions betw een the SDSS and Johnson photometric systems from 
ISmithet all (120021) . 

Table [3] lists the relevant results of mainly quasi- 
simultaneous optical/UV observations of 3C345, NRA0512 
and B3 1640+396 performed at Tzec Maun observatory and 
Swift UVOT. The Galactic extinction in the direction of these 
sources is between E(B -V) = 0.010 mag an d 0.014 mag accord- 
ing to the tables of .Schlegel et alj (Il998l) . Using the extinc- 
tion law by Carde lli et alJ (Il989l) a nd coefficients describing the 
UVOT filters (Romin g et alj|200^ . the following extinctions (in 
magnitudes) were obtained for the individual bands: 0.03-0.04 
(V), 0.04-0.06 (B), 0.05-0.07 (U), 0.07-0.09 (Wl), 0.09-0.13 
(M2) and 0.08-0.12 (W2). 

The source 3C 345 was monitored also at mm wavelengths 
(230 GHz) by the Sub-millimeter Array (SMA) on Mauna Kea 
for the GASR Additionally, a program conducted by Ann Wehrle 
(private communication) provided the data collected by the 
SMA for NRAO 512 and B3 1640+396. 

3. Results 

3.1. Localization of y -ray emission 

An unbinned spectral likelihood analysis was performed using 
integrated 20-month Ferm/-LAT monitoring data in the energy 
range of 0.1-lOOGeV. This energy range was chosen after a se- 
ries of different energy cuts were compared, to minimize the er- 
ror on the localization. In the IFGL catalog y-ray sources are 
typically localized with a 95% position eiTor radius of 0.02°- 
0.07° for sources with TS « 1000, but 0.06° -0.5° for TS -30 
dAbdo et al. 2010a). In cases of faint y-ray emitters it is difficult 
to obtain well-constrained localizations to make statistically sig- 
nificant associations, even more so for two faint y-ray emitters 
in close vicinity of each other 

A single point-source, 1FGLJ1642. 5+3947, is listed in the 
IFGL catalog around the analyzed region. Assuming a sin- 
gle point-source, the best-fit position using the 20-month LAT 
dataset is R.A. 16H2"24^ Dec. +39°48'27", with a 95% error 



Table 3. Quasi-simultaneous optical observations of 3C345, 
NRAO 5 12 and B3 1640+396 conducted by Swift UVOT and the 
200-mm telescope of Tzec Maun observatory. 



Source 



Obs. Date (*) 



Band Magnitude 



3C345: 



NRAO 512: 



B3 1640+396: 



2009-01-26 (4858) 
2009-02-13 (4876) 
2009-10-08(5113) 
2009-10-08(5113) 
2010-08-18 (5427) 
2009-10-08(5113) 
2010-08-06 (5415) 
2010-08-06 (5415) 
2009-01-26 (4858) 
2009-02-13 (4876) 



R 


16.48 + 0.02 


R 


16.92 + 0.03 


U 


16.44 + 0.06 


VWi 


16.28 + 0.03 


VWi 


16.47 + 0.02 


u 


18.61 + 0.04 


u 


17.71 +0.05 


VWi 


18.07 + 0.02 


R 


18.57 + 0.13 


R 


> 18.6+ 



Notes: Obs. Date (*) - Date of Swift-XRT observation 

(*=JD-2450000). 

Band - Optical filter used for photometric observations 

Magnitude ~ Observed optical magnitude in given filter band. The 

magnitudes are not corrected for Galactic extinction. 

t no detection 



circle radius of 0.037°. This procedure also provided the fol- 
lowing characteristic source parameters: photon indejqj F = 
2.44 ±0.03, integrated flux: (1.78 + 0.08)- lO^^phcm^^ s"', TS 
value of 2825. The estimated systematic uncertainty of the in- 
tegrated flux is 10% at 100 Me V, 5% at 500 Me V and 20% 
at lOGeV. After 20 months of data the localized position is 
shifted by 0.042° from its IFGL position toward NRAO 5 12, 
separated by 0.13° (east) from 3C345 and by 0.48°(west) from 
NRAO 5 12. This suggested a significant contribution of excess 
photons from the direction west of 3C 345. 

The unbinned spectral likelihood analysis was comple- 
mented by a pointlike analysis, similar to the method ap- 
plied to determine positions and error ellipses for the IFGL cat- 



* The photon spectral index F is defined as N(E) oc E '^, where N(E) 
is the y-ray photon flux as a function of energy E. 
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Table 4. The unbinned spectral likelihood results for the point- 
source localizations and spectra of 3C 345 and NRAO 512. 





3C345 


NRAO 5 12 




R.A. (J2000) 


16''43™0.24-^ 


16M0™44.4'' 




Dec. (J2000) 


+39°48'22.7" 


+39°46'12.0" 




Ar 


0.35' 


2.9' 




a (95%) 


0.0490° 


0.0735° 




b (95%) 


0.0442° 


0.0626° 


c 


O 


41.6 


-33.2 


■^ 


r(pi) 


2.49 ±0.02 


2.37 ±0.04 


c 


r(gt) 


2.45+0.05 


2.41 ±0.08 


s 


Flux (pi) 


1.22 + 0.10 


0.51 ±0.23 


Q 


Flux (gt) 


1.13±0.13 


0.67 ±0.13 




Pivot Energy 


949 


1598 




TS 


1076 


246 





Notes: 

R.A./Dec. ~ Right ascension/declination of the y-ray localization; 

Ar - distance to respective radio source position; 

a (95%) ~ major axis of 95% localization error ellipse; 

b (95%) - minor axis of 95% localization error ellipse; 

<t> - Rotation angle of the error ellipse East of North (degrees); 

r - spectral index of power-law fit, (pi) are the values obtained from 

the pointlike fit and (gt) are values from the gtlike fit; 

Flux - 20-month average flux (0. 1-100 GeV) determined from the 

power-law spectral fit in units of 10"' ph cm"^ s"' ; 

Pivot Energy - in MeV, this is the energy for which there is no 

correlation between the flux or normalization and spectral index 

uncertainties; 

TS - Likelihood test statistic value. 



al og. The app r oach o f pointlike is described in Section 4.2 
of lAbdo et al.l (l2010ah . This analysis revealed that a significant 
number of excess y-ray photons were detected from the vicinity 
of NRAO 5 12 after November 2009 (15 month). A fit with two 
point-source models in which source positions were not fixed, 
representing contributions by 3C345 and NRAO 5 12, yielded 
significant detections for the 20-month period (see Table |4|. 
Fig. [T] shows the radio positions of the three candidate sources 
together with the EGRET and IFGL localizations, as well as 
the improved localizations for the two y-ray sources found in 
the 20-month dataset. B3 1640+396 is 19.3' and 16.5' sepa- 
rated from the respective localizations of 3C 345 and NRAO 512. 
Both y-ray counterparts of 3C 345 and NRAO 512 now coincide 
with their respective radio positions (within the 95% error el- 
lipse). This suggests that both sources may be y-ray emitters. 
This is firmly established in Section [33] using the broadband 
temporal characteristics of the sources. There is no evidence 
for significant y-ray emission from the third candidate source, 
B3 1640+396. The results of pointlike were cross-checked 
with an iterative use of gtlike and gtfindsrc, which confirm 
the pointlike results within errors, see Table |4] A 2cr upper 
limit was calculated for B3 1640+396 by fixing each of the three 
candidate sources to their radio positions and assuming a F = 2.3 
spectral shape for B3 1640+396 and using the previously fitted 
F for 3C345 and NRAO 5 12. This resulted in a 20-month up- 
per limit for the integrated flux of 7.4-10"''phcm"^ s"'. Fig. |2] 
shows the observed binned spectra of the two sources together 
with their respective power-law fits. The pointlike best-fit pa- 
rameters are listed in Table |4] In Fig. [3] the respective 30-day 
binned y-ray light curves of the two localized sources are shown 
and are discussed below in more detail. 
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Fig. 1. y-ray counterpart localizations together with the posi- 
tions of the radio counterparts of the three candidate sources. 
Crosses mark the radio positions of 3C345, NRAO 5 12 and 
B3 1640+396. The large dotted orange circle denotes the 95% 
confidence eiTor l ocalization of EGRJ1642+3940 ( Casandiian 
& Grenier 2008). The dashed blue ellipse denotes t he 95% con- 
fidence error localization of IFGL J1642.5+3947 (l Abdo et al.] 
1201 Oah . The two solid green ellipses denote the 95% confidence 
error localization of 3C 345 and NRAO 512 presented in this pa- 
per. The confusing source Mrk501 is 2.5° east (left) and 4C 
+38.41 is 2.1° southwest of the field center. 



3.2. 20-month of y-ray monitoring 

3.2.1. Long-term variability of 3C 345 and NRAO 51 2 

To study the long-term variability, light curves were pro- 
duced for the y-ray counterparts associated with 3C 345 and 
NRAO 512 with 30-day binning. The 30-day binned light curve 
is shown in Fig. [3] For most of the data points NRAO 5 12 did 
not reach above a TS of 25, except for the last five months of the 
20-month period. This led to the conclusion that NRAO 5 12 did 
not produce a significant excess of y-ray emission over the local 
background during the first 15 months of LAT operation; thus it 
did not appear in the IFGL catalog. After 15 months the amount 
of excess photons from the vicinity of NRAO 5 12 slowly in- 
creased, making a detection over the integrated 20-month period 
possible. The average FermZ-LAT flux for the first 15-months of 
the monitoring period was (0.42 + 0.07) ■ 10"^ phcm"^ s"' com- 



,-2g-I 



for the last quarter of the 



pared to (1.1 +0.1)- lO^^phcm" 
monitoring period (months 15-20). 

The observed y-ray flux values reveal variability on 
timescales from days to months. To quantify this variability, the 
variability index { V) as defined through the;^'^ distribution (see 
lAbdo et al.ll2010ah was determined. It is computed from the 30- 
day integrated light curve data shown in Fig. [3] Whenever TS < 5 
the 2cr upper limit was computed and its error estimate for that 
interval was replaced with half the difference between that up- 
per limit and its value determined through the unbinned spectral 
likelihood analysis. In the absence of variability V is expected to 
follow a.x^ distribution with 19 (=Nint - 1) degrees of freedom. 
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Fig. 2. Twenty-month y-ray spectra of 3C345 and NRAO 5 12. 
The red curve represents the obtained respective average power- 
law fits. The red dots correspond to the pivot energies listed in 
Table H 



At the 99% confidence level the light curve is significantly dif- 
ferent from a flat one if V > 36.2. For 3C 345 a V of 1 15.3 was 
obtained, for NRAO 5 12 it was 37.7. Both sources are variable 
according to the 99% confidence interval with higher variability 
observed from 3C 345. 

In the shorter binned light curve of NRAO 512, ~63% of the 
derived values had a TS < 5 (i.e. non-detections) compared to 
-48% for 3C 345. Because NRAO 5 12 did not show a significant 
flux for most of the time on integrations of five days or less and 
with the observation that the flux of 3C 345 dominates, the model 
component of NRAO 512 was removed and light curves were 
produced with a single component fixed at the position of the 
y-ray locaUzation of 3C 345 listed in Table |4] This reduced the 
amount of noise in the light curve of 3C 345 and the number of 
values with a TS < 5 was only 14%. 

3.2.2. Short-term variability of 3C345 

For the remainder of this paper the discussion is based on light 
curves produced with a single point-source model component. 
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Fig. 3. y-ray Hght cui-ves of 3C 345 and NRAO 5 12 for the first 
20 months of Fenni-LAT monitoring, obtained through fitting 
two point-source model components with power-law spectra to 
the region of interest placed at the respective counterpart po- 
sitions. Values were determined from integrating observations 
over 30 days within an energy range of 0.1- 300 GeV. Blue 
filled circles are values obtained for 3C 345, red open circles are 
values for NRAO 5 12. Values with downward arrows represent 
2cr upper limits shown in place of detections with low signif- 
icance (TS<5). The dashed lines plot the respective 20-month 
average y-ray source flux. Errors are statistical only. 



The five-day binned light curve of 3C 345 is shown in Fig. |4]for 
comparison with observations at other wavelengths. 

To investigate shorter term variability in the observed y-ray 
emission of 3C 345, several prominent light-curve "flare events" 
were identified and parameters such as rise time, fall time and 
time of the peak were extracted using the five-day and two-day 
integrated light curves. These parameters can be used to con- 
strain the size of emission regions and cooling-times. Table |5] 
summarizes these results, listing characteristic parameters for 
all y-ray flare events for which the peak had a significance of 
TS > 25 in the two-day light curve. This corresponded to a flux 
greater than 3.0 ■ 10"^ ph cm"^ s ' . The time of onset and end of 
an individual flare event was determined through a change in the 
sign of the gradient at the rising and declining slope of the light 
curve around the peak of a flare event. 

Three dominant flare events of at least 5.0 • 10"^ ph cm"^ s"' 
and a significance of above 7cr were observed on two-day time 
scales, one during 2008 and two during 2009. They were labeled 
with I, II and III in Fig.|4] Several weaker flares were identified 
and labeled with numbers from 1 to 6. On average, the flares 
had a duration of 38 days a slightly faster rise than fall time and 
an average peak flux of (5.0 ± 0.5) • 10"^ phcm"^ s"' (two-day 
integrations centered on the peak). 

3.3. Identification of y-ray counterparts 

Following the y-ray outburst on October 2nd, 2009, Swift ob- 
servations were triggered on the candidate sources 3C345, 
NRAO 5 12 and B3 1640-1-396. X-ray spectra and optical mag- 
nitudes were obtained from these observations. A summary of 
relevant measurements is listed in Tables|2]and|3] 

The brightest source of the three candidates in X-ray as 
weU as in optical is 3C345. On January 26th, 2009, 3C345 
and B3 1640H-396 were detected in the optical R-band at re- 



Schinzel et al.: Identification of y-ray emission from 3C345 andNRA0 512 
Table 5. Characteristics of prominent y-ray events during the first 20 months of the Fermi-LAI observations. 



tpe-jk 


^stajt 


'^stop 


Af 


Afrise 


A?fall 




o m:ix 


(*) 


JD(t) 


JD(t) 


JD(t) 


days 


days 


days 




10-' p"; 




4734 


4709 


4754 


45 


25 


20 




5.5 ±1.5 


I 


4856 


4844 


4874 


30 


12 


18 




3.6 ±1.1 


1 


4916 


4904 


4944 


40 


12 


28 




4.0 ± 1.2 


2 


4978/4998 


4969 


5014 


45 


9 


16 


5.7 


± 1.3/5.4 ± 1.5 


II 


5020 


5014 


5044 


30 


6 


24 




4.5 ± 1.2 


3 


5056/5070 


5044 


5089 


45 


12 


19 


4.4 


± 1.4/4.5 ± 1.2 


4 


5104 


5099 


5119 


20 


5 


15 




9.1 ±1.6 


III 


5194 


5169 


5214 


45 


25 


20 




4.7 ±1.4 


5 


5278 


5264 






14 






3.6 ± 1.1 


6 


Average values: 






38 


13 


20 




5.0 ±0.5 





Notes: ?peak - in units of JD(t) = JD-2450000, pealc value was obtained from the two-day integrated light curve, all values have an error of ± 1 day; 



^ starts ^stop 
A?, Afriffi, Af| 



- in units of JD(t) = JD-2450000, all values have an error of ± 3 days; 



fall 



' in units of days, all values with error of ± 4 days; 



S"^" - two-day integrated flux value of the y-ray event peak value; 

(*) Labels are split into two sub-categories according to their peak flux and significance of detection, denoted by Roman and Arab numbers, see 



Fig-Hlfor details. 



spective magnitudes of 16.48 + 0.02 and 18.6 + 0.1. On October 
8th, 2009, observations of 3C345 and NRA0512 showed 
respective optical U-band magnitudes of 16.45 + 0.08 and 
18.64i0.04. Both NRA0512 and B3 1640+396 were a fac- 
tor of 7-8 fainter than 3C 345 at optical wavelengths. At X-ray 
energies, as seen from Tabled NRA0 512 and B3 1640+396 
had similar fluxes, which were a factor of 8-9 fainter than 
3C345. Archival Swift-UVOT and XRT data from observa- 
tions of NRA0512 made in January 2007 have an aver- 
age optical U magnitude of 18.0 ±0.1 and an X-ray flux of 
(8.3 + 2.3)-10-'^ergcm-^ s"', which are consistent with the re- 
cently observed values. Quasi-simultaneous mm observations 
between October and December 2009 determined the fluxes of 
3C345, NRA0512 and B3 1640+396 to be -3.0, -0.4 and 
~0.1-0.3Jy. The respective flux ratios are consistent with the 
ones in the optical/UV and X-rays. The mm radio flux density of 
3C 345 was at least a factor of 7 higher than for NRAO 512 and 
B3 1640+396. 

No significant excess of y-ray emission from NRAO 512 was 
observed during the first 15-months of the ferm/-LAT operation. 
At the end of 2009, NRAO 512 showed an increasing trend in its 
15 GHz VLBI rad io flux observed bv the MOJAVE team ( Lister 
et al. 2009n Indeed, the radio flux reached a low activity level 
around May 2009 and more than doubled until about July 2010 
with the mm-radio flux density reaching a peak value of over 
0.6 Jy in August 2010. A similar trend was observed in optical 
emission. Within one year, the optical flux increased by a factor 
of ~2.2 (see TableO. In X-rays a flux increase of at least 20-30% 
was observed, but because of the low statistics (less than 200 
photons) it is difficult to compare the observations. An increased 
number of y-ray photons was detected spatially consistent with 
NRAO 512 during the last quarter of the 20-month period. The 
excess is not high enough to affect the short-term bins of the 
light curve of 3C 345 significantly, however, as was shown in 
Section im the integrated 20-month dataset made it possible to 
detect a second point source west of 1FGLJ1642. 5+3947, which 
is consistent with the radio position of NRAO 512 (see Fig.[T]). 
Altogether this identifies NRAO 5 12 as a y-ray source. 



3.4. Multi-wavelength analysis of y-ray emission from 3C345 



|www.physics.purdue.edu/MOJAVE/sourcepages/1638+398.shtml| 



During 2009 nine optical observatories, listed in Section 12.31 
regularly monitored 3C 345, mainly as part of the GASP-WEBT 
program. This provided a densely sampled optical R-band light 
curve (Fig.|4]i. Visual comparison of the y-ray and optical light 
curves indicates a likely correspondence between at least four y- 
ray/optical flare events: Flare I, II, 4 and III. With the addition of 
X-ray and mm-radio monitoring data, the following aspects are 
of particular interest: 

1. A possible fast optical brightening of at least 0.3 mag 
with a duration of only a few days (Fig. |4]i was ob- 
served (JD 2454978), coincident with the peak of Flare II. 
Unfortunately, the optical monitoring had a gap between 
JD 2454953 and JD 2454977. At the same time the mm radio 
light curve showed a maximum of around 4.3 Jy. Moreover, 
a high X-ray flux was observed in the same period of the 
Flare II, followed by a decrease similar to what was observed 
in y and radio bands. Even if the maximum X-ray flux was 
observed two weeks before Flare II, the y-ray event lacks 
strictly simultaneous X-ray observations, thus we cannot ex- 
clude that a simultaneous X-ray flare was missed. 

2. The broad optical flaring episode observed around 
September/October 2009 had a total duration of about 
five months, with shorter timescale sub-structures that 
correspond to activity variations observed also at y- and 
X-rays. The peak in the optical was observed on September 
19th, 2009 (JD 2455094), and had no similarity in profile 
with events at other wavelengths (see Fig.|4]i. 

3. Ten days after the peak of the optical light curve, the strong 
y-ray Flare III was observed. The optical peak was fol- 
lowed by fast optical variability at a significant amplitude 
of 0.3-0.5 mag, with an elevated magnitude observed on 
September 29th on the same day as y-ray Flare III. 

4. Flares 1 and 3 can be related to the appearance of new 
features in the pc-scale radio jet of 3C 345 at 43 GHz (see 
ISchinzeletaljr2 010). 

5. All but Flares 5 and 6 occur before the significant excess of 
y-ray emission was observed near NRAO 512. However, for 
at least Flare 6, a correspondence in optical variability was 
found for 3C 345; unfortunately, a lack of multiwavelength 
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Fig. 4. Multiwavelength light curves of 3C 345 for the 20-month 
period discussed in this paper, from top to bottom for: a) y rays 
observed by Fermi LAT between 0. 1 and 300 GeV, values (filled 
circles) are five-day integrated, 2cr upper limits for five-day inte- 
grations are shown as light-blue crosses with downward arrows 
where TS < 5. Dark-blue filled circles have a TS > 25, the others 
are in the range 5 < TS < 25. In contrast to Fig. [3] the y-ray emis- 
sion was fitted by a single power-law component placed at the 
position of 3C345. b) X-rays observed by Swift XRT between 
0.3 and lOkeV. c) optical (R-band) observations performed by 
GASP and SAO. d) millimeter radio observations by the SMA 
at 230 GHz. The two orange downward arrows on top indicate 
the observation epoch for which a new feature was detected in 
the pc-scale radio iet of 3C 345, see Section 13.41 and Schinzel 
et al. (2010). Three dominant flare events of high significance 
observed on two-day time scales were labeled with I, II, and III. 
Several weaker events were identified and labeled with numbers 
from 1-6. 



data around Flare 5 makes it impossible to identify it with 
either 3C 345 or NRAO 512. 

To quantify the cross-correlation betw een light curves, a 
discre te correlation was calculated (DCF; Edelson & KrolikI 
11988 ^ using the 20-month Fenni-LAI light curve and the opti- 
cal R-band light curve between JD 2454680 (August 1st, 2008) 
and 2455299 (April 12th, 2010), as well as the radio SMA 
Ught curve between JD 2454721 (September 11th, 2008) and 
2455289 (April 2nd, 2010). Fig.|5](top) plots the resulting DCF 
comparing optical and y-ray light curves using a five-day y-ray 
integration and a ten-day timelag binning. This yields a number 
of distinct correlation peaks, of which the two most significant 
are discussed. In the first peak the optical leads the y-ray flux by 



(15 + 10) days with a correlation coefficient of 0.44 + 0.05. The 
second peak corresponds to the optical flux lagging the y-rays by 
110 + 5 days with a correlation coefficient of 0.77 + 0.06. To test 
the robustness of these peaks in the DCF, we artiffcially removed 
the points corresponding to the peak y-ray ffux from Flares II 
and III and recalculated the DCF. The same behavior was ob- 
served, with the y-ray ffux leading the optical by 20 days and 
lagging by 1 10 days. This shows that the results obtained do not 
solely depend on these two events. Finally, the DCF between the 
y-ray and the SMA light curve of 3C 345 was calculated, which 
yielded a noisy distribution of correlation peaks (also see Fig.|5]l. 
The strongest peak indicates that the radio emission leads the y 
rays by 120 + 5 days with a correlation coefficient of 0.83 + 0.30. 

The statistical significance of the cross-correlation peaks 
was investigated using Monte Carlo simulations, following 
Chatterjee et al. (2008) and Max-Moerbeck et al. (2010) and as- 
suming that the noise properties of the light curves can be de- 
scribed by a power-law power spectral density (PSD; esc l/f^^) 
with an exponent that depends on the particular energy band of 
the emission. 

In order to estimate the significance of the y-ray/optical cross 
correlation, an attempt was made to calculate and character- 
ize the PSD shape from the presented data directly, despite the 
limited length (20 months) and sampling (> 1 day) of the light 
curves. For the optical and y-ray data, the light curves available 
support a range of time scales corresponding to about one decade 
in the frequency domain. The respective range is even smaller 
for the more sparsely sampled radio data. The PSD s were cal- 
culated following Section 3. 1 in Uttle v et aP (l2002h . The result- 
ing power-law slopes, /?, are in the range of 1.2-0.7. Values of 
/? » 2.0 are not consisten t with the calculate d PSD shapes and 
can be ruled out. However. lUttlev et al.l (l2002h demonstrated that 
aliasing effects add to the power spectral slopes and could poten- 
tially ffatten the observed power spectrum. The relatively high 
noise and short time span of the data makes this range of yS rather 
unreliable. Monte Carlo simulations for these values of /3 yield 
extremely high cross-correlation signiffcances of >99.999% for 
the y-ray/optical light curves and >99.98% for the y-ray/radio 
case. It should be noted that a reliable determination of the PSD 
from observational data is usually done with a much larger dy- 
namic range (e.g. covering 4-5 orders of magnitude in th e fre- 
quency scale), which is also the case in lUttlev et al.l (l2002h . The 
radio, optical and y-ray data discussed in our paper do not pro- 
vide such an extended coverage in the frequency space. We adopt 
therefore a more cautious approach, exploring a range of power- 
law indices based on published values. 

In the y-ray regime ySy = 1.7 + 0.3 fo r BL Lacs and I3y = 
1.4 + 0.1 for ffat spectrum radio quasars (lAbdo et al.ll2010a) . 
Instead, for 3C 279, it has been found that 0rad w = 2.3 + 0.5 at 
14.5 G Hz and /^optical =" 1-7 + 0.3inflie R-ban d (IChatteriee et all 
12008 !). More recently IChatteriee efaLl (1201 Ih found an average 
optical power-law index of /^optical = 1 .6 ± 0.3 for a sample 
of six blazars in the R-band. Simulated light curves with PSD 
power-law indices in these ranges can provide reliable estimates 
on the signiffcance of the cross-correlations. The simulations use 
light curves of the same length as the data and simulated up to 
a maximum freq uency of 1 d ay'' using the method described 
by iTimmer & Koenid (Il995h . Each case uses 10000 pairs of 
light curves, which are sampled and cross-correlated in the same 
fashion as the observed light curves. From these simulated data 
the distribution of correlation coefficients for each time bin can 
be obtained for the case of non correlated time series. Fig. |5] 
presents the 3cr (99.7%) signiffcance levels for all tested PSD 
shape combinations. 
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The y-ray/optical correlation has a significance of higher 
than 3cr (99.7%) for the peak at a timelag of 1 10 days in most 
of the tested cases. It is the most significant peak in the DCF 
and it has a moderately high confidence level (> 98.5%) even 
for the least significant case {Py - 2.5 andjSopticai = 2.5; highest 
dashed line in the top panel of Fig. |5]l among the fi values con- 
sidered. The y-ray /optical peak at a timelag of about -20 days 
has a much lower confidence level (> 67.6%) in the least sig- 
nificant case and it remains below the 3cr level at all times. The 
y-ray/radio cross-correlation at timelag -120 days has a confi- 
dence level of > 90.7% in the least significant case {fiy - 2.5 and 
Aadio = 3.0). 

The source 3C 345 is by far the brightest of the three quasars 
in the y-ray region studied during the presented time period. 
In all bands it was consistently brighter than NRA0 512 and 
B3 1640+396. The improved Fermi position combined with 
multi-wavelength activity and a moderate significance of cor- 
related optical/y-ray variability, using a plausible range of fi in 
the y-ray and optical PSDs, strongly suggest an identification of 
3C 345 as y-ray source and the main contributor to the observed 
y-ray flux in IFGL J1642.5-H3947. 



4. Discussion and summary 

Properties of the y-ray emittin g region reported in the first 
three-month Bright Source List jAbdo et al.ll2009a|) as well as 
in th e first 11-month Fenni-LAT source catalog (lAbdo et alj 
1201 O a.b) were compared to the 20-month dataset presented here. 
Localizations of the sources in this region were improved over 
the previously reported values, leading to the detection of two y- 
ray point sources, one consistent with IFGL J1642. 5+3947 and 
a new y-ray source ~ 0.4° west of it. 3C 345 and NRAO 512 were 
identified as the counterparts for these y-ray sources. Their y-ray 
spectra have a spectral slope of F = 2.4-2.5 typical for that of flat 
spectrum radio quasars (.Abdo et al]|2010bh . No significant spec- 
tral break was observed i n the spectrum of either source, unhke 
in the case of 3C 454.3 dAckermann et aljl2010h . however, the 
statistics at the highest energies are poor, and the presence of a 
spectral break at E>10GeV is not ruled out. 

The y-ray emission from 3C345 was identified based 
on an improved y-ray counterpart localization and multi- 
wavelength activity including correlations of moderate signif- 
icance found between the optical and radio variability, us- 
ing a plausible range of /3 in the y-ray and optical PSDs, 
and major y-ray events observed by Fermi LAT during its 
first 20 months of operation. EGRET observed 3C 345 be- 
tween 1991 and 1994 during a time of enhanced radio activ- 
ity. An upper limit of 2.5 ■ 10"^ phcm"^ s"' for th e viewing pe- 
riod S eptember 1 2th - 19th, 1991, was reporte d (iFichtel et alJ 
[Illi). However, I Casandiian & Grenied ( |2008|) report the de- 
tection of a y-ray source that could be consistent with 3C 345 
in the period April 23rd - May 7th, 1996 with a y-ray flux of 
(3.5 + 0.8)- lO^^phcm^^s"', but it was noted that the associ- 
ation with 3C 345 remains unclear because of a possible con- 
fusion with Mrk501. In the Fermi-era, with an increased sen- 
sitivity and the possibility of long-term integrations, 3C 345 
was already detected during the first three months of opera- 
tion at a level of ~ 18cr and after 20 months it reached over 
30cr. However, owing to the unclear association of 3C 345 with 
EGR J1642+3940, the three-month data were falsely associated 
with CLASS J1641+3035 (B3 1640+396), which, as shown in 
this paper, did not produce a significant excess of y-ray emission 
above the detection threshold of Fermi-hP^. Thus the simplest 
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Fig. 5. Correlation between the five-day integrated y-ray light 
curve and the optical (top) and radio (bottom) light curves of 
3C345. The correlation is obtained by calculating the discrete 
correlation function (DCF) with ten-day binning. A positive 
time-lag corresponds to the y-ray variations preceding the vari- 
ations in the other two bands. The significance of the correla- 
tion is illustrated by 3cr (99.7%) significance contours calculated 
from simulated, uncorrected data with different slopes, /?, of the 
power spectral density (PSD). The contours are drawn for the 
PSD slopes p - 1.0, 1.5, 2.0, and 2.5 (indicated by progres- 
sion of colors, from light gray to black). For the simulated data, 
values of the correlation coefficient exceeding +1 are an inher- 
ent feature of the method, resulting from dif ferences in overlap 
of the time series at different time lags (cf. lEdelson & KrolikI 
19881; iTimmer & Koenig 1995, and the description of the ap- 
plied method in Section lT^l l. 



explanation for the previous non-detection of 3C 345 is that its 
emission was just below or at the detection limit of EGRET. 

These findin gs end the decade-lon g debate whether 3C 345 
is y-ray loud (lUnwin et al.l Il997l) . The identification of 
1FGLJ1642.5+3947 with 3C345 is further supported by the 
connection of radio-em ission of the parsec-s cale jet with the ob- 
served y-ray emission (iSchinzel et al.ll2010l) . Here not only in- 
dividual flares were correlated with radio events, but an unam- 
bigious long-term trend in y-ray emission was observed, which 
matched the trend of the radio flux density light curve of the in- 
ner jet at 43 GHz. The trend and correlated radio events relate to 
a distance of ~10pc in the source frame and questions existing 
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scenarios that place the y-ray emission at sub -pc distances from 
the central engine (e.g., Tavecchio et al. TOW). This supports the 
arguments given by Unwin et al. ( 1997) that the observed y-ray 
emission is related to the relativistic jet on pc scales. 

The cross-correlation found between y-ray and optical light 
curves at 1 10 days had a moderate confidence level of > 98.5%, 
whereas the lag at -20 days had a confidence level of only 
> 67.6% using a plausible range of /3 in the y-ray and optical 
PSDs. If typical values for the power spectral density (PSD) 
shapes are assumed (J3y = 1.5 and ySopticai = 1-5; second to 
lowest dashed line in the top panel of Fig. [S]), the confidence 
level remains low at > 89.3%. The weak correlation found be- 
tween the y-ray and 230 GHz radio light curves had a worst-case 
confidence level of >89.3%. Assuming a possibly more typi- 
cal combination of PSD shapes (J3y = 1.5 and ySi-adio - 2.0), a 
confidence level of > 96. 1 % is obtained. If this weak correlation 
holds, it could suggest a connection between the radio emission 
and production of y-ray s. A timelag of ~-120 days suggests the 
mm-radio emission precedes the r-ravs. |Pushkarev et al.l (1201 Ol) 
showed for a larger sample of radio-loud quasars that on av- 
erage the 15 GHz radio emission lags ~1.2 months behind y- 
rays. The result presented here indicates that between 15 GHz 
and 230 GHz it is possible to d irectly observe the y -ray emitting 
region, as already implied by lUnwin et al.l (Il997l) . which sup- 
ports the argument of y-ray production in the relativistic jet on 
pc scales. However, in contrast to these, a more confident cross- 
correlation result of y-rays leading by ~80 and ~50 days was re- 
cently obtained for tw o prominent BL La c type objects, OJ28 7 
(lAgudoetal.ll2011ah and AO 0235-1-164 ( Agudo et al.ll2011bh . 
respectively. Doubts on the indicated radio/y-ray correlation of 
3C 345 remain in this problem and should be investigated more 
thoroughly with a longer and denser sampled lightcurve at mm 
wavelengths. 

The y-ray emission from NRAO 5 12, a high-redshift quasar, 
was identified based on a correlated increase in multiwavelength 
emission from radio up to GeV energies between 2009 and 2010. 
At monthly timescales the source started to be significant around 
November 2009. This increase of y-ray flux led to the first detec- 
tion of the now identified new y-ray counterpart west of 3C 345. 
The LAT flux quoted in this paper for the first 15 months was 
below the EGRET detection threshold, offering a plausible ex- 
planation of why EGRET did not see emission from NRAO 512, 
except possibly during a y-ray flare event. 

The radio, optical and X-ray flux ratios between 3C 345, 
NRAO 512 and B3 1640-1-396 could be used to gauge their indi- 
vidual y-ray photon contributions and might be able to provide 
an indication of false associations for y-ray photons from this 
part of the sky. Multiwavelength monitoring, in addition to the 
continuous all-sky monitoring by Fermi LAT, will help to im- 
prove the localization of y-ray emission from this region and is 
necessary for a deeper variability study for any of these sources 
in the future. Ultimately this might lead to the detection of y-ray 
emission from B3 1640+396 as well. 
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